numerical solution of partial differential equation